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ABSTRACT 

Examination of the process of silicon burning, the burning stage that leads 
to the production of the iron peak nuclei, reveals that the nuclear evolution 
is dominated by large groups of nuclei in mutual equilibrium. These quasi- 
equilibrium (QSE) groups form well in advance of the global Nuclear Statistical 
Equilibrium (NSE). We present an improved "minimal" nuclear network, 
which takes advantage of quasi-equilibrium in order to calculate the nuclear 
evolution and energy generation while further reducing the computational 
cost compared to a conventional a-chain network. During silicon burning, 
the resultant QSE-reduced a network is twice as fast as the full a network it 
replaces and requires the tracking of only half as many abundance variables, 
without significant loss of accuracy. When the QSE-reduced a network is used 
in combination with a conventional a network stretching from He to Si, the 
combined a7 network provides an accurate approximation for all of the burning 
stages from He burning to NSE, while tracking only 7 abundances. These 
reductions in computational cost and the number of species evolved make the a7 
network well suited for inclusion within hydrodynamic simulations, particularly 
those in multi-dimension. 
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Subject headings: methods: numerical — nuclear reactions, nucleosynthesis, 
abundances — stars: evolution — supernovae: general 

1. Introduction 

Nucleosynthesis calculations perform two functions in astrophysical models. The most 
important function, in the context of the evolution of the encompassing hydrodynamic model, 
is the calculation of the rate of thermonuclear energy generation. Where nucleosynthesis 
occurs, it is often the dominant source of energy, making careful determination of this 
rate of energy generation vital to accurate hydrodynamic modeling. The second function 
of a nucleosynthesis calculation is to determine the evolution of the nuclear species. 
For nuclear astrophysics and our understanding of the origin of the elements, this is of 
paramount importance. Often the best comparisons of model calculations with observations 
also depend on detailed knowledge of the abundances of nuclear species. Unfortunately, 
accurate calculation of the nuclear evolution is computationally expensive. For this reason, 
it is common to include within a hydrodynamic model only a limited approximation 
for nuclear burning which estimates the rate of energy generation. The detailed nuclear 
evolution is then calculated at a later time using the thermodynamic trajectories from the 
hydrodynamic model. Such a scheme is referred to as post-processing nucleosynthesis. The 
success of such a post-processing scheme depends largely on the accuracy with which the 
limited approximation which is included within the dynamical calculation can estimate 
the rate of energy generation. Frequently this inline nuclear burning calculation is a small 
network, trading a limited ability to follow the nuclear evolution for speed. 

Tracking the nuclear evolution from helium burning through to NSE requires a network 
that includes nuclei from a-particles to Zn. Silicon burning presents a particular problem 
as material proceeds from silicon to the iron peak not via heavy ion captures but through 
a chain of photodisintegrations and light particle captures. The minimal nuclear set which 
can follow this evolution is the set of a-particle nuclei; a, ^^C, ^^O, ^°Ne, ^^Mg, ^^Si, ^^S, 
36 Ar, ^°Ca, ^^Ti, ^^Cr, ^^Pe, ^^Ni, ^^Zn. For convenience we will label this full set and 
refer to its abundance as Y-^. Silicon burning in fact presents a larger problem, as the 
nuclear flow from silicon to the iron peak nuclei does not generally procede through nuclei 
with N=Z, especially when significant neutronization has occurred ( [Hix fc Thielemann 1996 
: henceforth HT96). In some models, however, such compromise is made necessary by the 
computational limitations. 
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The nature of nuclear network calculations has been extensively reviewed elsewhere 
(see, e.g., [Woosley 1986| , [I'hielemann et al. 1994| or [Arnett 1996| ), so we will discuss it only 



briefly here. From a set of nuclear abundances, and the reaction rates for the reactions 
which link them, one can calculate the time derivatives of the abundances, Y. With these 
derivatives, the abundances of the included nuclei, at time t, and over timestep At, are 
evolved according to the finite difference prescription 

Y{t + At) - Y{t) _ Y{t) (explicit) 

Y{t + At) [implicit) . 

For the stiff set of non-linear differential equations which form most nuclear networks, a 
fully implicit treatment is most successful. Solving the implicit version of Eq. ^ is equivalent 
to finding the zeros of the set of equations 

^> + Ai).I<i±M^_y(i + A„.0. (2) 



This is done using the Newton- Raphson method (see, e.g., [Press, et al., 1992]), which is 



based on the Taylor series expansion of Z{t + At), with the trial change in abundances 
given by 

AY=(m±^Yz, (3) 

where dZ/dY is the Jacobian of Z. Iteration continues until Y{t + At) converges. 

The conclusion of silicon burning with nuclear statistical equilibrium (NSE) offers 
a hint to a more efficient means of evolving nuclear abundances, since the equilibrium 
distribution reduces the number of independant variables which must be tracked. NSE 
grows from several groups of nuclei which form local equilibria. The existence of these 
quasi-equilibrium groups frees one from accounting for the changes in the abundances of 
each of the members of these groups, as the changes within the groups can be accounted 
for by the changes of a few crucial abundances. Previous authors, most notably [Bodansky^ 



Clayton fc ii'bwler (1968 : henceforth BCF), have attempted to use quasi-equilibrium to 
simplify the calculation of silicon burning. Having noted the computational difficulties 
encountered by P?ruran et al. (1966)| in their pioneering nuclear network study of silicon 



burning, BCF sought to model silicon burning with a single quasi-equilibrium group, 
stretching from ^^Mg and ^'^Si through the iron peak nuclei. BCF postulated that the 
downward fiow from this group occurred via ^^Mg(7, a)^''Ne and that the fiow of a-particles 
downward through ^°Ne, ^^O, and ^^C conserved the abundances of these nuclei. From 
these assumptions, they derived expressions for the abundances of these light a-particle 
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nuclei and for the downward flux of a-particles which depended only on the abundances 
of a-particles and ^'^Si. Though this model was moderately successful in replicating the 
abundance pattern found in the solar system from A=30 to A=60, it had clear limitations 
when compared to the network calculations of TCG, particularly at early times when the 
assumption of conservative flow was most suspect. [Woosley, Arnett, fc Clayton (1973 : 



henceforth WAC) concluded that such efforts were at best marginally successful, in part 
because as matter cools the photodisintegration reactions freeze out well before the particle 
captures cease, thereby destroying quasi-equilibrium. WAC found a further complication 
in that, at early times, their models exhibited two separate QSE groups, one around 
silicon and the other containing the iron peak nuclei. WAC did demonstrate that, with 
the assumption of QSE, reasonable estimates of the time required to merge the two QSE 
groups and the time to reach silicon exhaustion could be calculated. HT96 demonstrated 
additional complications as both the path of the nuclear flow from the silicon group to the 
iron peak group and the membership of these groups depended strongly on the amount 
of neutronization. HT96 further demonstrated that for material which has undergone 
significant neutronization, the separation between the silicon and iron peak groups persists 
through a significant part of the nuclear evolution. In spite of these complications, pix 

Thielemann (1998 : henceforth HT98) showed that QSE could provide a reasonable 
estimate of the abundances of many species during sihcon burning, even under explosive 
conditions, where the photodisintegration reactions freeze out before their corresponding 
particle captures. Arnett 1996| provides a more complete history of the use of QSE. 



2. The QSE-reduced Network 

The objective of the QSE-reduced a network is to evolve F-^ during silicon burning, 
and calculate the resulting energy generation, in a more efficient way. Under conditions 
where QSE applies, the existence of the silicon and iron peak QSE groups allows calculation 
of the abundances of these 14 nuclei from a reduced set of 7 abundances. For the members 
of the silicon group (^^Si, ^^S, ^^Ar, ^°Ca, ^^Ti), the individual abundances can be calculated 
from the assumption of QSE by 

CA^Z^ A-28 

YqseA^Z) = , (4) 

where we have defined 
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for later convenience and 



The functions G{^Z) and B{^Z) are the partition function and binding energy of the 
nucleus ^Z, Na is Avagadro's number, fc^ is Boltzmann's constant, and p and T are the 
density and temperature of the plasma. is the abundance of free a-particles and the 
integer [A — 28)/4 is the number of a-particles needed to construct "^Z from ^'^Si. Similarly, 
the abundances of the members of the iron peak group (^^Cr, ^^Fe, ^^Ni, ^°Zn) can be 
calculated by 

YqseA^Z) = ^|^^(''Ni)y„^ . (6) 

Thus, under conditions where QSE applies, Y-^ can be expressed as a function of the 
abundances of the reduced nuclear set TZ, defined as [a, ^^C, ^^O, ^°Ne, ^^Mg, ^^Si, ^^Ni], 
and it is therefore sufficient to evolve only Y^. It should be noted that WAC and HT96 
have shown that ^^Mg is ordinarily a member of the silicon QSE group during silicon 
burning, leaving only 6 independant abundances, but for ease in integration with the 
conventional nuclear network which we will discuss in §4, we will evolve the ^^Mg abundance 
independant ly (see [Arnett 1996| ). 



While Y^ is a convenient set of abundances for the calculation of Y-^, it is not the 
most efficient set to evolve, primarily because of the complexity of calculating the necessary 
time derivatives and their Jacobian. It is much easier to calculate the time derivatives of 
each QSE group as a whole. To this end, we define 9^= [F„g, Y{^^C), Y{^^0), Y(20Ne), 
Y(24Mg), YsiG, Yf,g], where 

Y.a = Y^+ E ^^y^+ E ^^Y., 

iGSi group iGFe group 

YsrG = E (7) 

i£Si group 

YpeG = E 

i&Fe group 

Physically, Ysic and Ypec represent the total abundances of the silicon and iron peak 
QSE groups, while Y^g represents the sum of the abundances of free a-particles and those 
a-particles required to build the members of the QSE groups from ^^Si or ^^Ni. 

Employing the group abundance set Q also reduces the number of reactions whose flux 
must be calculated since quasi-equilibrium allows one to ignore the reactions among the 
members of the QSE groups. Unfortunately, the rates for this reduced set of reactions are 
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still functions of Y-^ . For example, the nuclear flux into the iron peak group, 

YpeG = p{cxv)Y^Y C^Ti) - XYC'Ct) , (8) 

depends not on the group abundances, Yac, Ysic and Yp^c, but on the abundances of free 
a-particles, ^^Ti, and '^^Cr. Since these time derivatives are not easily expressed in terms of 
one must calculate Y^ from Y-^ or alternately, using Eqs. ^ and from Y^. Applying 
Eqs. H and |^ to Eq. ^ results in 

y.ec = P{-')^v(-si)yi - A§p^y(»N,)y;- . p) 

Thus, to evolve Y^ over time, one must repeatedly solve for Y'^. By applying Eqs. ^ and ^ 
to Eq. 1^ one may write Y^ as a function of Y''^, with seven independant equations for the 
seven abundances of Y^ implying a unique relation between Y^ and Y^. Fortunately, with 
the exception of Y^, the relations between Y^ and Y''^ are linear. 

Further complicating the calculation, however, is the need in Eq. |^ for the Jacobian of 
Z. This requires knowledge of dY^ / dY^ , which can not be calculated directly since Y^ can 
not be expressed in terms of Y^ . Instead, we have been successful using the chain rule, 

dY^ dfs dY^ , , 

= ^ 10 

gys dY'^dY^ 

to calculate the Jacobian. Analytically, the first term on the righthand side of Eq. |10| is 
easily calculated from the sums of reaction terms. For example, from Eq. ^ we can see that 
the non-zero terms of dYpeG/ dY''^ are 



9F(28Si) ' 'C(28Si) " 9F(56Ni) C(56Ni) 

Calculation of the second term on the righthand side of Eq. |T0| is more complicated, 
requiring implicit differentiation of the definition of Y^ with respect to F^, using Eq. |^. 
As an example, in order to calculate dYpeGl dYp^G, we differentiate Eq. |^ with respect 
to YpeGi resulting in 3 equations for three unknowns, dYa/ dYpeGi dY i^^'^i) / dYpeGi and 
dY i^^Wi) / dYpeGi while the other 4 terms of dY^ / dYpec are manifestly zero. Solving these 
3 equations results in the 3 terms 
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dY.^o M%,M%!^''' ^^^^ 



9y(56Ni) F(56Ni) / M\' M 



1 + 



where we have used the following definitions for simplification 

ieSi Group \ ^ / jgPe Group V ^ / 



- ,13) 



Multiplying each term in Eq. |T2| by the respective component of Eq. |TT] and summing 
produces dYp^cl dYp^G- 



3. Silicon burning with the QSE-reduced a network 

In this section we will demonstrate the accuracy with which the QSE-reduced a 
network duplicates the results of the full 14 element a network for silicon burning. Our 
first examples are nucleosynthesis calculations occuring under constant temperature and 
density. While these calculations provide the least challenging comparison, they also allow 
comparison with NSE, which should represent the final abundances of these calculations. 
Figure |I| offers comparison of the mass fractions of the 7 independant species; a-particles, 
^^C, ^^O, 20Ne, 24Mg, and the sihcon and iron peak groups, as evolved by the QSE-reduced 
and conventional a networks for silicon burning at 5 x 10^ K and a density of lO^gcm"^. 
Apart from an early enhancement by the QSE-reduced network of the iron peak mass 
fraction (20% after 10~^ seconds), these mass fractions typically agree to within 1%. Since 
the nuclear energy release depends linearly on the abundance changes, differences in small 
abundances have little effect on the nuclear energy store. In this case, the difference in the 
rate of energy generation calculated by the two networks is < 1% at 10~^ and 10~^ seconds. 
This difference is significantly smaller than the variation, shown by both networks, in the 
rate of energy generation between timesteps, with e typically declining by 5% per timestep 
over this interval. 

With such good agreement between the respective abundances of free a-particles and 
the two QSE group abundances, significant variations in abundance among the individual 
members of the QSE groups can only result from deviations from QSE. At early times, 
the small abundances within the iron peak reduce the accuracy of QSE at predicting the 
individual abundances of members of the iron peak group. Much of the enhanced mass 
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Fig. 1. — Evolution of the independant miclcar mass fractions for constant thermodynamic 
conditions, Tg=5 and p = lO^gcm"^. The sohd hnes display the evolution due to a 
conventional a network, the circles show the evolution by the QSE-reduced a network. 
The Silicon group mass fraction is the sum of the mass fractions of ^^Si, ^^S, '^^Ar, ^°Ca, and 
44rp| rjnj^^ j^^^ group mass fraction is the sum of the mass fractions of ^^Cr, ^^Fe, ^^Ni, and 
60Zn 
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fraction of the iron peak nuclei at early times, seen in Fig. |l|, is due to the QSE-reduced 
network's emphasis on heavier nuclei at the expense of ^^Cr. After an elapsed time of 
10~^ seconds, the average mass of the iron peak nuclei, Apec, is 49.2 according to the 
conventional network and 52.6 according to the QSE-reduced network. As a result, the 
abundances of "^^Cr and ^^Fe calculated by the QSE-reduced network are 38% and 164% 
of their conventional network values, while ^^Ni and ^^Zn are 16 times more abundant 
than the conventional network predicts. As the iron peak nuclei become more abundant, 
QSE provides a better estimate of the relative abundances within the group, reducing 
such discrepancies. By the time the iron peak nuclei represent a significant portion of the 
mass, the differences in the abundance predictions for all nuclei are only a few per cent. 
Comparison of columns 2 and 3 of Table |1] show that after an elapsed time of 10~^ seconds, 
^^Cr displays the largest difference, with the QSE-reduced network's abundance equal 
to 97% of the abundance calculated by the conventional network. In the context of the 
network evolution, this small difference is less than half of the abundance change for ^®Cr 
that occurs over the following timestep. As each network reaches its respective equilibria 
the abundance predictions of these networks continue to differ by at most 3%, even among 
the nuclei with the smallest abundances. Not surprisingly, in view of these small abundance 
differences, the difference in the total energy released by these networks is less than .1%. 
Comparison of the network abundances with abundances calculated from NSE show a 
similarly low level of difference. It should be noted that these NSE abundances, calculated 
assuming that the a nuclei are the only available states in the equilibrium, differ noticably 
from NSE calculations which include all nuclei, even for = 0.5. 

For higher temperatures, large abundances of free a-particles reduce the importance of 
the silicon and iron peak nuclei. As Column 4 of Table |^ reveals, even the NSE abundance 
distribution which represents the end point of silicon burning may be dominated by free 
a-particles. In spite of the larger role played by free a-particles, the QSE-reduced a network 
successfully replicates the nuclear evolution during silicon burning, because the nuclei of 
the silicon and iron peak groups still obey QSE. Fig. ^ compares the nuclear evolution of 
the independant mass fractions from the QSE-reduced a network to their conventional 
counterparts, with Tg=6 and p = lO^gcm^^. Here too, the iron group mass fraction is 
initially over-predicted by the QSE-reduced network, though to a lesser extent (7% after 
10~^ seconds). This over-prediction also subsides as the iron peak abundances increase, 
and hence their adherence to QSE improves. As in the previous case, despite these small 
differences in abundance, the agreement between the rates of energy generation is better 
than 1%. Examination of Table ^ reveals that the equilibrium abundances for the individual 
nuclei calculated by the QSE-reduced network also agree well with those calculated by the 
conventional network and by NSE. Even for the smallest abundances, the differences are 
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Table 1: Comparison between the QSE-reduced and conventional network abundances after 
1.06 xlO~^ seconds with Tg=5 and p = lO^gcm"^ 



Nucleus 


net 




Y 

qse 




^He 


1.67x10" 


-5 


1.68x10" 


-5 




2.59x10" 


-6 


2.58x10" 


-6 


160 


4.09x10" 


-6 


4.09x10" 


-6 


20Ne 


3.29x10" 


-7 


3.28x10" 


-7 




4.11x10" 


-5 


4.11x10" 


-9 


28Si 


1.88x10" 


-2 


1.89x10" 


-2 


32S 


8.40x10" 


-3 


8.36x10" 


-3 


36 Ar 


2.07x10" 


-3 


2.04x10" 


-3 


^"Ca 


1.27x10" 


-3 


1.24x10" 


-3 




1.52x10" 


-5 


1.48x10" 


-6 


48cr 


5.63x10" 


-5 


5.45x10" 


-5 


52pg 


2.67x10" 


-4 


2.69x10" 


-4 


56Ni 


1.11x10" 


-3 


1.13x10" 


-3 


60Zn 


6.65x10" 


-8 


6.73x10" 


-8 



Table 2: Comparison of network abundances at feezeout with NSE for Tg—6 and p — 
10''gcm~3 



Nucleus 














^He 


1.23x10" 


-1 


1.24x10" 


-1 


1.25x10" 


-1 


12^ 


1.03x10" 


-6 


1.04x10" 


-6 


1.04x10" 


-6 


16Q 


1.88x10" 


-6 


1.91x10" 


-6 


1.89x10" 


-6 


20Ne 


3.54x10" 


-8 


3.62x10" 


-8 


3.56x10" 


-8 




4.08x10" 


-6 


4.18x10" 


-6 


4.08x10" 


-6 


28si 


1.28x10" 


-3 


1.32x10" 


-3 


1.28x10" 


-3 


32S 


1.21x10" 


-3 


1.23x10" 


-3 


1.20x10" 


-3 


36 Ar 


7.03x10" 


-4 


7.13x10" 


-4 


7.00x10" 


-4 


^•^Ca 


7.60x10" 


-4 


7.64x10" 


-4 


7.55x10" 


-4 


44rjJ 


4.16x10" 


-5 


4.14x10" 


-5 


4.12x10" 


-5 


48cr 


2.70x10" 


-4 


2.71x10" 


-4 


2.67x10" 


-4 


52pg 


1.43x10" 


-3 


1.42x10" 


-3 


1.41x10" 


-3 


56Ni 


5.13x10" 


-3 


5.06x10" 


-3 


5.06x10" 


-3 


60Zn 


2.98x10" 


-6 


2.92x10" 


-6 


2.94x10" 


-6 



- 11 - 




-10 -8 -6 -4 -2 

Log (Time) s 



Fig. 2. — Evolution of the indcpcndant nuclear mass fractions for constant thermodynamic 
conditions, Tq=6 and p = lO^gcm"^. The solid lines display the evolution due to a 
conventional a network, the circles show the evolution by the QSE-reduced a network. 
The Silicon group mass fraction is the sum of the mass fractions of ^^Si, ^^S, '^^Ar, ^°Ca, and 
44rp| rjnj^^ j^^^ group mass fraction is the sum of the mass fractions of ^^Cr, ^^Fe, ^^Ni, and 
60Zn 
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less than 



70. 



While examples of silicon burning calculations under constant conditions are instructive, 
if the QSE-reduced a network is to replace the conventional a network it must be shown 
to be accurate under changing thermodynamic conditions. Of particular importance is the 
ability to model explosive silicon burning. As a representative analytic model for silicon 
burning occurring as a result of shock heating, we will consider a mass zone which is 
instantaneously heated by a passing shock to some peak initial temperature Tgj and density 
Pi and then expands and cools adiabatically. Following the approximation introduced 
by [I'bwler fc Hoyle (1964)| , the hydrodynamic expansion timescale (equal to the free fall 



timescale) is 

thd = (247rG'p)"^/^ = 446p6 ^^Vs, (14) 

with the time dependence of the density and temperature of this radiation dominated gas 
given by 

p(t) = piexp (-t/THD) , 

T^{t) = T9,exp(-t/3rHZ5), (15) 

assuming the adiabatic exponent 7 = 4/3. For initial conditions we chose conditions 
identical to the examples discussed above. 

Figure ^ shows the nuclear evolution for an example of this explosive burning model 
with Tgj = 5 and pi = lO^gcm"^. Over the first millisecond, the evolution portrayed 
here is virtually identical to that of Figure |I|; however by the time one millisecond has 
elapsed, the temperature has dropped to 4.9 xlO^ K, slowing the reactions which are 
turning silicon into iron peak elements. This cooling, which drops Tg below 4 after 9 ms 
and below 3 after 22 ms, freezes out the nuclear reactions before NSE is reached, resulting 
in incomplete silicon burning, as discussed by WAG. In this case, the freezeout leaves nearly 
equal amounts of silicon group and iron peak group elements. Throughout most of the 
evolution in this example, the agreement between the mass fractions as evolved by the 
QSE-reduced a network with those evolved by its conventional counterpart is comparable 
to that demonstrated under constant thermodynamic conditions. Columns 2 and 3 of 
Table |^ compare the abundances after 9 ms have elapsed, with Tg nearing 4. In this case, 
as for that displayed in Table |I|, the largest relative error (5%) is in the abundance of 
^^Cr. These small differences in abundance result in small differences in the accumulated 
nuclear energy release, approximately .5% to this point. By this time, adiabatic cooling has 
greatly reduced the rate of energy generation from its peak of more than 10^^ ergs g^^ s^^ 
to roughly 10^^ ergs g^^ s~^. Though the absolute difference in the rate of energy generation 
as calculated by the 2 networks has declined from ~ 10^^ to 10^^ ergsg"^ s~^, the relative 
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Fig. 3. — Evolution of the independant nuclear mass fractions under adiabatic expansion 
with Tgi—5 and pi — lO^gcm"^. The sohd hues display the evolution due to a conventional 
a network, the circles show the evolution by the QSE-reduced a network. The Sihcon group 
mass fraction is the sum of the mass fractions of ^^Si, ^^S, ^^Ar, ^°Ca, and ^^Ti. The Iron 
group mass fraction is the sum of the mass fractions of ^^Cr, ^^Fe, ^^Ni, and ^°Zn 
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difference has grown to 10% as Tg nears 4. Fortunately this difference is neghgible, since 
nuclear energy release has virtually ceased, with < .2% of the total energy release remaining. 

Though energetically unimportant, nuclear reactions continue, resulting in significant 
changes in the smaller abundances. As pointed out by WAC, and discussed by HT98 
in the context of a large nuclear network, the continued cooling results in the gradual 
freezeout of the photodisintegrations responsible for QSE. Typically, a photodisintegration 
only represents a significant flux if its Q value is less than 30 ksT ~ 2.6 Tg MeV. The a 
captures which link silicon to nickel typically have Q values of 7-8 MeV, thus QSE begins 
to fail as Tg approaches 3. However this same decline in temperature also reduces the rate 
of charged particle capture reactions, greatly reducing the amount of nucleosynthesis which 
occurs after Tg drops below ~ 3.5. The group abundances of the silicon and iron peak 
groups (which account for 99.9% of the mass), as calculated by the QSE-reduced a network 
after 18 ms have elapsed (Tg=3.3), differ by less than 1% from those of the conventional 
a network at the same point in time. Comparison of Columns 4 and 5 of Table ^ reveals 
larger variations among individual abundances, most notably, significant under estimation 
by the QSE-reduced network of the ^*Cr and ^^Fe abundances, with a small, compensatory 
over estimate of the ^^Ni abundance. These variations, factors of 2 and 4 for ^^Fe and ^^Cr, 
respectively, and 3% for ^^Ni signal the breakdown of the iron peak QSE group. With the 
steep decline in temperature and density, the flux upward from ^^Fe in the conventional 
network is no longer sufficient to provide the reduction in abundance which QSE and the 
sharply declining abundance of free a-particles requires. 

As the temperature and density continue to decline, so to does the free a-particle 
abundance. Column 6 of Table ^ details the abundances after 255 ms have elapsed, with Tg 
having dropped to .01, and all abundances having frozen out. Comparison of colums 4 and 
6 reveals that the decline of the free a-particle abundance by 10~^ is the largest abundance 
variation beyond 18 ms. Since the more abundant species have effectively frozen out by 
the time Tg approaches 3.5, comparison of columns 5 and 6 reveal that the predictions of 
the QSE-reduced network, frozen near Tg=3.5, also provide good abundance estimates, in 
addition to the excellent estimates of the rate of energy generation discusssed earlier, in 
spite of the freezeout. 

Figure ^ shows the evolution of another example of this explosive silicon burning 
model with Tgj = 6 and pi = lO^gcm"^. The lower density in this case results in a slower 
expansion, with Tg reaching 5, 4, and 3 after 77, 172, and 293 milliseconds, respectively. 
The nuclear evolution is similar to that portrayed in Figure ^ for the first 10 milliseconds, 
after which the declining temperature favors the iron peak group at the expense of lighter 
nuclei. Though the rate of cooling is relatively slow, the temperature drops too quickly for 
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Table 3: Comparison of network abundances for T9j=5.0 and pi= 10^ gem 



Time (ms) 




8.77 






17.7 




255 




^9 




4.07 






3.29 




0.01 




Nucleus 










^' net 












^He 


7.90x10" 


-7 


7.82x10" 


-7 


1.04x10" 


-8 


1.01x10" 


-8 


1.94x10" 


-14 




1.96x10" 


-7 


1.96x10" 


-7 


3.99x10" 


-8 


3.23x10" 


-8 


3.90x10" 


-8 


160 


7.34x10- 


-7 


7.39x10" 


-7 


5.26x10" 


-7 


5.07x10" 


-7 


5.27x10" 


-7 


20Ne 
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Fig. 4. — Evolution of the independant nuclear mass fractions under adiabatic expansion 
with Tgi—6 and pi — lO^gcm"^. The sohd hues display the evolution due to a conventional 
a network, the circles show the evolution by the QSE-reduced a network. The Sihcon group 
mass fraction is the sum of the mass fractions of ^^Si, ^^S, ^^Ar, ^°Ca, and ^^Ti. The Iron 
group mass fraction is the sum of the mass fractions of ^^Cr, ^^Fe, ^^Ni, and ^°Zn 
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the large a-particle abundance to be completely incorporated into heavier nuclei, resulting 
in a-rich freezeout (see |Arnett, Truran, fc Woosley 1971| , and [Woosley, Arnett, fc Clayton 



(1973| ). Even with this large overabundance of a-particles, the QSE-reduced network tracks 



the group abundances and rate of energy generation as reliably as it did under constant 
conditions until Tg approaches 3.5. By the time 223 ms have elapsed (T9=3.5), the largest 
relative difference, occuring in the small (~ 10~^) abundance of the silicon group, has grown 
to 9%. With only such small errors in abundance, it is not surprising that the difference in 
the rate of energy generation remains less than 1%. There is greater disparity between the 
individual abundances as calculated by the QSE-reduced network and those calculated by 
the conventional network. 

As discussed by HT98 for a more complete nuclear network, in cases of a-rich freezeout 
the large a-particle abundance results in a relatively large nuclear flow upward into the 
silicon group. This flow gradually disrupts the QSE groups by increasing the abundance 
of the less massive members of the group while their QSE abundances drop sharply. 
Fortunately for our use of QSE to reduce the network, because of the large free a-particle 
abundance, these less massive members have only small abundances and thus the effects 
of this disruption are small. Comparison of abundances calculated by the conventional a 
network after 170 ms have elapsed (Tg ~ 4) with those of its QSE-reduced counterpart 
(columns 2 and 3 of Table ^) shows the beginning of this process, with the abundances 
of ^^Si and ^^S much larger than QSE would predict. In spite of this breakdown in QSE, 
the abundances predicted for the more abundant members of the group and, as a result, 
the rate of energy generation, agree well. However as the temperature continues to drop, 
the disruption of the QSE groups by the large free a-particle abundance grows. By the 
time Tg drops to 3 (columns 4 & 5 of Table |[), the under prediction by QSE of the smaller 
group abundances also effects the iron peak group. The abundances of the dominant nuclei, 
a, ^^Ni and ^°Zn, and hence the energy production, are still in good agreement, with the 
energy generation predictions differing by less than 1%. 

Of particular importance, the abundance of ^°Zn is within 2% of that predicted by QSE. 
With a Q value of 2.7 MeV, ^^Ni -f a <-> ^°Zn + 7 remains balanced until significantly later 
than the reaction pairs connecting silicon to nickel. As the remaining photodisintegrations 
freeze out, the continued capture of the large abundance of a-particles results in a large 
number of a captures and significant abundances changes, the most prominent of which 
is the conversion of ^^Ni to ®°Zn. As late as an elapsed time of .35 seconds, with Tq=2.6 
and y(^°Zn) = 1.5 x 10~^, the abundance calculated for ^°Zn by the QSE-reduced network 
is within 5% of that predicted by the conventional network. As a result of this continued 
validity of QSE for the abundance of ^°Zn, the energy generation rate predicted by the 
QSE-reduced a network at an elapsed time of .35 seconds is only 11% greater than that 
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Table 4: Comparison of network abundances near freezeout for T9i=6.0 and pj=10 gem 



Time (s) 




.174 






.290 




2.67 




T9 




3.98 






3.02 




0.01 
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^' iicL 
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-2 
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-2 
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-2 
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-2 
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-9 
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-9 
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-6 
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-10 


2.46x10" 
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1.90x10" 


-8 


6.44x10" 


-9 
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-6 
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4.15x10- 


-5 
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calculated by the full a network. As the temperature continues to drop, the continued 
a captures disrupt this last balanced reaction pair, converting even more ^^Ni into ^°Zn. 
Column 6 of Table ^ shows the abundances after 2.7 seconds have elapsed, by which time 
all of the abundances are frozen. While virtually all of the abundances have grown at 
the expense of those of free a-particles and ^^Ni, the most significant is a doubling of the 
abundance of ^°Zn after .35 seconds. Though the QSE-reduced network can not track the 
nuclear evolution to this point, the relatively small amount of energy released, less than 1% 
of that released over the first .35 seconds, makes this only a minor weakness. 

4. The (x1 network 

While the preceding sections demonstrate the success of the QSE-reduced a network 
as a replacement for a conventional a network during silicon burning, it can not accurately 
calculate the evolution of Y'^ when QSE does not apply. Prior to silicon burning, the 
primary abundances of interest are those of a, ^^C, ^^O, ^°Ne, ^^Mg, and ^^Si. This subset 
of nuclei, which we will call 7?.', is identical to 7?., less ^^Ni. Thus evolution of a conventional 
a network for the 7 nuclei of the set 71. can model these earlier burning stages nearly as well 
as the larger a network which evolves the entire set JF, with the only difficulty being small 
abundances of nuclei heavier than silicon. In combination with the QSE-reduced a network, 
one could hope to globally replace the full a network with a smaller network which evolves 

using QSE where it applies. The principle uncertainty in the use of this combined 
network, which we have dubbed a7, is the the nature and timing of the transition to QSE. 

The first necessary condition for QSE is that the temperature be in excess of 3 x 10^ K, 
because, as discussed in §3, for lower temperatures the rates of photodisintegration are much 
smaller than the corresponding captures. An additional requirement for QSE to provide 
a good estimate of the abundances of heavier elements is that a significant fraction of the 
matter be composed of nuclei with A > 24. As long as the material is principally composed 
of lighter nuclei, like ^^C, ^^O, or ^°Ne, the production and destruction of free a-particles is 
dominated by reactions among these nuclei, rendering suspect QSE abundance predictions 
for the silicon and iron peak group nuclei based on these free a-particle abundances. Once 
these conditions are met, the abundances of the set TZ' must be mapped in a suitable way 
into the initial abundances for the set Q which the QSE-reduced a network evolves. 

One approach to this transition seeks to conserve the individual abundances of the 
sets TZ and TZ' across the transition, setting Y'^' = Y^, and then calculating the initial 
values for Y^ according to Eq. |^. Figure |^ displays the energy generation under constant 
conditions (79=3 and p = lO^gcm"^) for the a7 network using this individual abundance 
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Fig. 5. — Comparison of the Energy generation by the full a network under constant 
thermodyamic conditions Tg^S and p — lO^gcm"^ to the aJ network, assuming the 
individual abundance transition. The solid curve is the result of the conventional 14 element 
a network. The dotted curve shows the results of the al network, with the transtion to QSE 
turned off. The short dashed, long dashed and dot dashed curves display the results of the 
al network, with the transition to QSE occuring when approximately 30%, 60% and 90% of 
the mass is in nuclei heavier than A — 24, respectively. 
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Log Time (s) 

Fig. 6. — Comparison of the evolution of the independant mass fractions by the full a 
network under constant thermodyamic conditions T9=3 and p = lO^gcm"^ to the a7 
network, assuming the individual abundance transition. The solid curves are the result 
of the conventional 14 element a network. The dotted curves show the results of the a7 
network with the transtion to QSE turned off. The short dashed, long dashed and dot 
dashed curves display the results of the a7 network, with the transition to QSE occuring 
when approximately 30%, 60% and 90% of the mass is in nuclei heavier than A — 24, 
respectively. 
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(lA) approach to the transition. Figure ^ shows the group abundances for this same lA 
case. For Figures 5-8, the sohd curves are the result of the conventional 14 element a 
network and the dotted curves show the results of the 7 element a network alone. The 
short dashed, long dashed and dot dashed curves display the results of the combined a7 
network, with the transition to QSE occuring when approximately 30%, 60% and 90% of 
the mass is in nuclei heavier than A = 24, respectively. The spikes in the short dashed, 
long dashed and dot dashed curves seen in Fig. ^ result from errors in mass conservation 
at the transition. Treating Y'^' as results in an over estimate of the silicon and iron 
peak group mass fractions at the expense of ^^O, producing the spike in energy generation. 
This overabundance of iron also causes the excess energy production around 100 s and the 
corresponding dearth around 10^ s. 

If one examines the a-particle abundance (shown in Fig. ^) in the vicinity of 10~^ s, one 
sees that the 7 element a network, and its a7 network descendants, begin to overproduce 
free a-particles. In the conventional 14 element a network, these a-particles are capturing 
on ^^Si, populating the Si group at the expense of ^'^Si, an option not allowed the smaller 
network. It is these overlarge abundaces of ^'^Si and free a-particles which cause the errors 
in mass conservation when QSE is first applied. A successful transition requires awareness 
that the abundance of ^^Si in the set TZ' represents the total of all abundances above ^'^Mg. 
In terms of Y^ , the capture of an a-particle on ^^Si does not change Yac or Ysic- If we 
consider Y^ when nuclei heavier than ^*Si are not evolved, the sums in Eq. |^ disappear, 
leaving Y^g = Ya, YsiG = Y{^^Si), and Yp^c = ^^(^^Ni) = constant. Thus Y''^' can also 
be mapped directly to Y^ . Since and Y^ have the same normalization rule, mass is 
conserved. This approach seeks to conserve the abundances of the groups across the QSE 
transition, hence we refer to this as the group abundance (GA) transition. Figure ^ shows 
the energy generation, again for burning at a constant Tg=3 and p = lO^gcm"^, using this 
approach to transition. The results agree more closely with the conventional a network 
for all three values of the transition mass fraction, with the best results occuring for a 
transition mass fraction around 50%. The mass fractions predicted by this approach, shown 
in Figure ^, are also more reliable, though there are still significant discrepancies among the 
smaller mass fractions, particularly near the transition to QSE and for a transition mass 
fraction of 90%. 

5. Conclusion 

In this paper, we have demonstrated the ability of QSE to reduce the number of 
independant nuclear abundances which need to be evolved to track a-chain nucleosynthesis 
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Fig. 7. — Comparison of the Energy generation by the full a network under constant 
thermodyamic conditions Tg^S and p — lO^gcm"^ to the aJ network, assuming the group 
abundance transition. The solid curve is the result of the conventional 14 element a network. 
The dotted curve shows the results of the al network, with the transtion to QSE turned off. 
The short dashed, long dashed and dot dashed curves display the results of the aJ network, 
with the transition to QSE occuring when approximately 30%, 60% and 90% of the mass is 
in nuclei heavier than A — 24, respectively. 
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Fig. 8. — Comparison of the evolution of the independant mass fractions by the full a 
network under constant thermodyamic conditions Tg—3 and p — lO^gcm"^ to the a7 
network, assuming the group abundance transition. The solid curves are the result of the 
conventional 14 element a network. The dotted curves show the results of the a7 network 
with the transtion to QSE turned off. The short dashed, long dashed and dot dashed 
curves display the results of the a7 network, with the transition to QSE occuring when 
approximately 30%, 60% and 90% of the mass is in nuclei heavier than A — 24, respectively. 
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for silicon burning. An a network reduced in this fashion is capable of providing reasonable 
estimates of the nuclear energy generation and elemental production, even in cases where 
thermodynamic variations result in incomplete silicon burning or a freezeout rich in free 
a-particles. Computational, this reduction in the number of independent nuclear abundance 
variables offers two advantages. First, within the nuclear network calculation, the smaller 
number of variables results in a smaller matrix to build and solve. This speeds the 
nucleosynthesis calculation by a factor of 2 in the case of the a network, with the potential 
for much greater speed increases for larger networks. Second, the reduction in the number 
of nuclear variables also reduces the number of equations which must be solved within 
the encompassing hydrodynamic model, a matter of particular concern for models with 
large numbers of zones or grid points, especially multi-dimensional studies. The relative 
importance of these 2 advantages depends on the size of the network. As a rule of thumb, 
a PPM hydrodynamic code evolving 45 nuclear abundances spends roughly half its time 
on the network solution ( |Arnett 1998| ). Thus for small networks, like the a network, the 



reduction in necessary number of hydrodynamic equations is the paramount advantage. 

In an effort to streamline the computation necessary for advanced nuclear burning 
stages, we have combined a smaller a network (for calculating the nuclear evolution from 
helium burning through oxygen burning) with this QSE based method. The resulting 
combined network, which we have dubbed a7, can be used in place of the full 14 element a 
network from helium burning through to NSE without significant errors in energy generation 
or nucleosynthesis. Such a combination offers the same reductions in the computational 
cost discussed above for silicon burning for all burning stages beyond H burning. In a future 
paper [Hix, Freiburghaus, fc Thielemann (1999)| , we will apply the method of QSE-reduction 
to larger nuclear networks, where the potential for improvement in speed and size is even 
greater. 
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